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ABSTRACT 

A  boundary  value  problem  for  the  Tricomi  equa- 
tion is  studied  in  connection  with  transonic  gas  dynamics. 
The  transformed  equation 

Au  +  3Y  "y  "  ° 

in  canonical  coordinates  is  considered  in  the  complex  domain 
of  two  independent  complex  variables.   A  boundary  value 
problem  is  then  set  by  prescribing  the  real  part  of  the 
solution  on  the  boundary  of  the  real  unit  circle. 

The  Dirichlet  problem  in  the  upper  unit  semicircle  with 
vanishing  values  of  the  solution  at  Y  =  0  is  solved  explicitly 
in  terms  of  the  hypergeometric  function  for  the  more  general 
Euler-Poisson-Darboux  equation.   An  explicit  representation 
of  the  solution  is  also  given  for  a  mixed  Dirichlet  and 
Neumann  problem  for  the  same  equation  and  domain. 

Reflection  rules  are  given  for  these  solutions  which 
permit  one  to  extend  them  from  the  upper  to  the  lower  unit 
semicircle.   The  transonic  boundary  value  problem  is  solved 
by  expressing  the  solution  as  a  linear  combination  of  these 
two  types  of  solutions. 


1.    INTRODUCTION 

The  Tricomi  equation 


(1.1)  yu    +  u 

■^  XX     yy 


is  an  example  of  a  partial  differential  equation  of  mixed 
type,  that  is,  an  equation  which  is  partly  elliptic  and 
partly  hyperbolic.  This  equation  has  been  widely  studied 
because  of  its  relation  to  transonic  gas  dynamics   [10] . 

The  Tricomi  problem  is  a  model  of  a  well  posed  problem 
for  equation  (1.1).   We  recall  [10]  that  it  is  formulated 
in  a  domain  which  is  intersected  by  the  parabolic  line 
y  =  0.   In  the  elliptic   half -plane  the  domain  is  bounded 
by  a  simple  curve  ending  at  two  points  on  the  parabolic 
line.   In  the  hyperbolic  half-plane  the  domain  is  bounded 
by  two  intersecting  characteristics  that  start  at  these 
two  points  of  the  parabolic  line.   The  Tricomi  problem  is 
well  posed  because  the  boundary  values  of  the  solution  are 
prescribed  on  the  whole  elliptic  boundary  but  only  on 
one  of  these  characteristics. 

An  important  problem  in  transonic  gas  dynamics  is  the 
computation  of  the  flow  past  a  given  airfoil  with  speed  at 
infinity  close  to  but  smaller  than  the  speed  of  sound.  In 
this  situation  the  flow  over  the  upper  side  of  the  airfoil 
may  become  supersonic.  The  partial  differential  equations 
of  motion  are  then  of  mixed  elliptic  and  hyperbolic  type. 


The  problem  of  calculating  continuous  transonic  flow  past 
such  an  airfoil  was  shown  by  C.  Morawetz  [20]  to  be  not 
well  posed. 

Motivated  by  the  method  of  complex  characteristics,  a 
transonic  boundary  value  problem  was  developed  in  [9]  for  a 
certain  domain  obtained  by  analytic  extension  to  complex 
values  of  the  independent  variables.   Complex  conjugate 
characteristic   coordinates  C  and  n  were  introduced.  In 
these  coordinates  the  stream  function  \p{^,E,)    is  real  in  the 
subsonic  domain  but  complex  otherwise.   The  new  boundary 
value  problem  is  set  by  prescribing  the  real  part  of  \p{E,,I,) 
on  the  boundary  of  the  unit  circle.   This  problem  has,  as 
yet,  not  been  shown  to  be  well  posed,  although  its  numerical 
solution  was  found  successfully.   Our  contribution  will  be 
to  show  that  a  corresponding  problem  for  the  Tricomi  equation 
is  well  posed. 

Under  the  change  of  coordinates 

9  „2     3 

4  Y   =  y   ,    X  =  X 

equation  (1.1)  can  be  written  in  the  canonical  form 

(1.2)  Au  +  -^   u^  =  0  . 

We  consider  this  equation  in  the  four-dimensional  space 
where  the  variables  X  and  Y  are  complex.   Points  in  the 
upper  unit  semicircle  of  the  real  (X,Y) -plane  correspond 
to  points  in  the  elliptic  half-plane  of  the  variables  x  and  y, 


which  we  shall  call  the  Tricomi  variables.   Points  in  the 

lower  unit  semicircle  of  the  real  (X,Y) -plane  correspond 

to  a  complex  surface  in  the  space  of  the  x  and  y  coordinates. 

A  solution  u(X,Y)  of  (1.2),  real  in  the  upper  unit 
semicircle,  is  continued  analytically  to  the  lower  unit 
semicircle  over  a  path  which  avoids  the  singular  locus  Y  =  0. 
The  boundary  value  problem  in  which  we  are  interested  con- 
sists in  prescribing  the  real  part  of  u(X,Y)  on  the  boundary 
of  the  real  unit  circle  in  the  (X,Y) -plane.   A  solution  of 
this  boundary  value  problem  will  be  given  that  can  be 
described  formally  as  follows: 

Two  types  of  solutions  for  equation  (1.2)  can  be  found. 
One  of  them  is   of  the  form 

(1.3)  u^(X,Y)  =  Y^/^  f(X,Y^)  , 

where  f(X,Y  )  is  an  analytic  function  even  in  Y.   The  other 
is  of  the  form 

(1.4)  U2(X,Y)  =  g(X,Y^)  , 

2 
where  g(X,Y  )  is  again  analytic  and  even  in  Y. 

The  Dirichlet  problem  in  the  upper  unit  semicircle  with 

vanishing  boundary  values  at  Y  =  0  can  be  solved  by  a  solution 

of  the  type  (1.3).   The  mixed  Dirichlet-Neumann  boundary 

value  problem  in  the  same  domain,  with  vanishing  normal 

derivative  on  Y  =  0  and  prescribed  values  of  the  solution  on 

the  remaining  boundary,  can  be  solved  by  a  solution  of  the 


type  (1.4).  By  analytic  continuation,  these  solutions  can 
be  extended  to  the  lower  unit  semicircle.  The  solution  u„ 
will  still  be  real.   It  will  be  of  the  form 

(1.5)        u„(X,-Y)  =  g(X,Y^)  ,    Y  >  0  . 


(1.6)   u^(X,-Y)  =  |y|^/^  e2TTi/3  f(x,Y^)  ,    Y 


>  0 


when  a  proper  branch  is  chosen. 

In  Chapter  5   it  will  be  shown  that  the  solution  u  of 
our  boundary  value  problem  can  be  expressed  in  a  unique 
way  as 

(1.7)         u(X,Y)  =  u^(X,Y)  +  U2(X,Y)  , 

where  u,  and  u„  are  of  the  form   (1.3)  and  (1.4)  respectively. 
A  uniquely  invertible  relationship  between  the  boundary  values 
of  u  and  those  of  u,  and  u„   is  obtained  from  the  evenness 
of  f  and  g   as  functions  of  Y.   Using  (1.5)  and  (1.6)  we 
can  write  this  relationship  as 

u(X,Y)  =  |y|^/^  f(X,Y^)  +  g(X,Y^)  ,    Y  >  0  , 
and 
Re{u(X,Y)}  =  -  I  |y|^/-^  f(X,Y^)  +  g(X,Y^)  ,    Y  <  0  . 

This  shows  that  the  boundary  values  of  u,  and  u„  on  the 
upper  unit  semicircle  can  be  chosen  so  that  the  real  part 
of  u  satisfies  the  prescribed  boundary  conditions  on  the 
whole  perimeter  of  the  unit  circle. 
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2.    AN  APPLICATION 

As  we  have  mentioned  in  Chapter  1,  our  boundary  value 
problem  has  been  motivated  by  the  problem  of  calculating 
shockless  transonic  airfoils  by  the  method  of  complex  char- 
acteristics described  in  [9] .   In  this  chapter  we  give  a 
description  of  the  method,  specifically  of  those  aspects 
which  motivate  our  problem. 

For  an  understanding  of  the  problem  it  is  important  to 
find  a  suitable  set  of  coordinates  in  which  it  can  be  stated. 
We  start  by  expressing  the  equations  of  motion  in  the  systems 
of  coordinates  that  are  relevant.   The  boundary  value  problem 
of  reference  [9]  is  then  stated.   Afterwards  we  give  a  des- 
cription of  transonic  paths  of  integration  that  are  required. 
We  consider  the  geometry  primarily  in  the  case  of  the  Tricomi 
equation  and  the  Tricomi  gas.   This  is  a  hypothetical  gas 
whose  equation  of  state  is  such  that,  in  pertinent  coordi- 
nates, the  equation  for  the  stream  function  becomes  the 
Tricomi  equation  (see  [4]). 

2.1.  Equations  of  Motion  and  the  Hodograph  Transformation 

Let  us  consider  the  equations  of  plane  flow  for  an 
inviscid  compressible  fluid.   We  assume  steady,  isentropic 
and  irrotational  flow.   Consider  an  ideal  polytropic  gas 
whose  equation  of  state  is  given  by  p  =  p ''^ ,  so  the  process 
is  isentropic.   The  variables  p  and  p  are  respectively  the 
pressure  and  density,  and  the  constant  y  is  the  adiabatic 


exponent  of  the  gas. 

For  steady,  irrotational  flow  the  Bernoulli  law  becomes 

In  this  equation  q  is  the  speed,  c  is  the  speed  of  sound, 
and  c^  is  the  critical  speed,  which  is  constant.   Let  u  and  v 
be  the  components  of  the  velocity  vector  in  the  (x,y) -plane. 
The  continuity  equation 

(2.2)  (pu)^  +  (pv)   =  0 
and  the  irrotationality  condition 

(2.3)  u   -  V  =  0 

y    X 

permit  us  to  introduce  a  stream  function  'j^  and  a  velocity 
potential  ^P .      Equations  (2.2)  and  (2.3)  can  then  be  trans- 
formed, using  the  Bernoulli  law  (2.1),  into  second  order 
partial  differential  equations  for  (t>  or  i(j  .  We  have,  for 
instance,  the  equation 

(2.4)  (c^-u^)(})    -  2uv(})    +  (c^-v^)4)    =  0  . 

^xx       ^xy  ^yy 

The  equation  (2.4)  is  quasilinear.   It  can  be  reduced 
to  a  linear  equation  by  use  of  the  hodograph  transformation. 
The  hodograph  transformation  consists  of  introducing  u  and  v 
as  the  independent  variables.  In  the  hodograph  plane  we  use 
polar  coordinates  defined  by 


The  equation  (2.4)  is  transformed   into  the  linear  system 

(2.5) 

*e  =   p  ^q    ' 

where  M  =  q/c  is  the  local  Mach  number. 

Equations  (2.5)  are  known  in  the  literature   as  the 
Chaplygin  equations.   They  are  of  the  hyperbolic  type  for 
supersonic  flow,  M   >  1,  and  of  the  elliptic  type  for 
subsonic  flow,  M   <  1.   The  ordinary  differential  equation 
for  the  characteristics  of  (2.5)  is 

(1-M^)  (dq)^  +  q^(d8)^  =  0  . 


Introducing  characteristic  coordinates  E,,t],      we  obtain  the 
canonical  system 


(2.6) 


1^  -  /i^  q^ 


qe^  +  A^   q^  =  0 


in  the  hyperbolic  domain. 

We  can  extend  the  concept  of  characteristics  to  subsonic 
flow  by  considering  analytic  extension  of  the  function  ip   into 
the  four-dimensional  domain  of  the  complex  variables  q  and 
e,  or  x  and  y  [10].   The  solutions  of  (2.6)  are  then  the 
complex  characteristics  and  the  characteristic  coordinates 
^  and  n  become  independent  complex  variables.   In  terms  of 
the  complex  characteristic  coordinates  we  obtain  the  canonical 


equations 


(2.7) 


h  =  ^Vc  ' 


\  =  i^-^ 


for  <p   and  ip  where 


/^^ 


An  analytic  function  of  a  characteristic  coordinate  is 
a  new  characteristic  coordinate.   Thus,  there  is  a  conformal 
mapping  free  that  can  be  used  to  transform  the  characteristic 
coordinates  conveniently.   Introducing  the  modified  speed  h, 
defined  by 

(2.8)  ^=  A^  ^    , 

we  can  write  the  equations  of  the  characteristics  as 

log  h  +  10  =  n / 
(2.9) 

log  h  -  ie  =  ^. 

By  use  of  a  conformal  transformation  new  characteristic 
coordinates   s  =  s(C)  and  t  =  t(n)  can  be  defined  as 

(2.10)       s  =  h(q)e"^^  ,    t  =  h(q)e^®  . 

Then  a  further  mapping  s=e  ^   ,t=e  ^      can  be  intro- 
duced to  transform  the  flow  onto  the  unit  circle. 
The  system  (2.7)  can  be  transformed  putting 

q   


(2.11)     X  =  9  ,    Y  =  log  h  =     /1-M' 


into  the  linear  equation  in  canonical  form 

(2.12)  *XX  +  *YY  -^  ?  S  *Y  =  " 

for  ip ,    with  T  =  T  .   After  another  change  of  coordinates 
defined  by 


da  =  p   g 


the  Chaplygin  system  transforms  into  the  equation 


(2.12)  K{a)^^ 

with 

K(a)  =  T   =  — 2~   • 

P 

This  may  be  viewed  as  a  generalized  Tricomi  equation. 

The  characteristic  coordinates  s  and  t  introduced  in 
(2.10)  have  the  property  that,  for  M  <  1,  s  =  t  if  and 
only  if  q  and  6  are  real.   Two  complex  coordinates  with  this 
property  are  called  conjugate  coordinates  [1] .   More 
generally  let  x  and  y  be  any  two  complex  variables,  and  let 

* 
z=x+iy,     z   =x-iy. 

We  see  that  z   =  z  if  and  only  if  x  and  y  are  real.  There- 
fore we  call  z  and  z*  conjugate  coordinates. 

Starting  with  the  complex  coordinates  x  and  y  there  is 
a  simple  way  of  generating  new  conjugate  coordinates.   Let 
f (z)  be  an  analytic  function  which  is  real  on  the  real  axis. 
Define  ^  =  f(z)  and  n  =  f(z*).   Then  C  and  n  are   conjugate 
coordinates.   For  if  x  and  y  are  real,  then  z*  =  z  ,  so  by 


the  Schwarz   reflection  principle   f(z  )  =  f(z)  =  f(z) 
and  ^  =  n.   The  converse  is  also  true,  provided  that  f  is 
one  to  one  in  the  domain  in  question. 

2.2.  The  Inverse  Problem  for  Supercritical  Airfoil  Design 

The  use  of  the  hodograph  transformation  is  appropriate 
for  solving  inverse  problems  in  transonic  gas  dynamics.  By 
the  inverse  method  a  continuous  flow  can  be  computed  in  terms 
of  hodograph  coordinates.  Consider  the  streamlines  in  the 
physical  plane.   If  the  streamline  4;  =  0  is  closed  and  not 
self  intersecting, it  can  be  taken  as  a  profile  which 
generates  the  flow. 

A  new  idea  is  given  in  [9]  to  solve  the  problem  of  air- 
foil design.   This  idea  is  to  map  the  unknown  domain  of  the 
flow  in  the  hodograph  plane  conformally  onto  the  unit  circle 
of  the  plane  of  one  characteristic  coordinate.   The  speed  q 
on  the  airfoil  is  prescribed  as  a  function  q(s)  of  the  arc 
length  s.   What  is  intended  is  to  have  the  equations  of  motion 
solved  for  the  stream  function  in  such  a  way  that  the  boundary 
of  the  unit  circle,  for  subsonic  points,  is  a  streamline  and 
furthermore  its  image  in  the  hodograph  plane  acquires  the 
prescribed  values  of  q  over  the  profile. 

More  general  conjugate  characteristic  coordinates  E,    and 
n  can  be  introduced  by  means  of  the  formulas 

(2.13)         s  =  e'^^^     ,      t  =  e^^   , 

where  the  map  function  f (C)  is  an  analytic  function  defined 
10 


in  the  unit  circle  \E,\    <    1,  and  s  and  t  have  been  defined  in 
(2.10).   The  characteristic  coordinates  E,    and  n  are  conjugate 
in  the  subsonic  domain,  so  points  in  the  real  subsonic  domain 
correspond  to  points  ^   =    r\.      Beyond  the  sonic  line  points 
C  =  1  do  not,  however,  correspond  to  real  supersonic  flow. 
The  stream  function  4;{C,f)  is  thus  real  in  the  subsonic 
domain  but  complex  otherwise. 

A  new  boundary  value  problem  was  formulated  in  [9]  for 
the  stream  function  ij;(^,ri).   Let  ip  be  a  solution  of  the 
canonical  system  (2.7)  that  is  allowed  to  have  a  pole  together 
with  a  logarithmic  singularity  in  the  hodograph  plane.  The 
boundary  condition 


Re 


|4^(^,e)}  =  0  ,      1^1=1 


is  imposed.   For  this  purpose  a  suitable  branch  of  the 
analytic  function  i>{^,r\)    has  to  be  specified. 

Paths  which  are  complex  conjugates  of  each  other  are  used 
to  compute  the  solution  in  the  real  subsonic  region  [1,9,14]. 
Given  the   two  initial  characteristic  planes   ^  =  Cq   and 
T]   =   T]      ,    a   path  5  =  ^(a)  is  defined  in  the  plane  n  =  tIq  • 
The  conjugate  path  n(T)  is  defined  in  the  characteristic  plane 
C  =  C_  as  the  set  of  points  n  =  E,  (a)  .   Consider  the  rectangle 
formed  by  taking  the  Cartesian  product  of  these  paths.  The 
diagonal  of  the  rectangle  is  given  by  the  set  of  points 
(C/H)  =  (C/i)  and  they    correspond  to  real  subsonic  points 
until  the  sonic  line  is  reached. 
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To  compute  points  below  the  sonic  line  special  paths 
have  to  be  defined  that  avoid  the  complex  sonic  surface  M  =1, 
where  the  canonical  equations  become  singular.   To  achieve 
this  the  paths  have  to  be  separated  so  that  they  are  no 
longer  conjugate.   An  analytic  continuation  of  the  solution 
around  the  sonic  line  that  is  appropriate  for  the  boundary 
value  problem  is  defined  by  such  integration  paths. 

The  geometry  in  the  space  of  the  canonical  coordinates 
X  and  Y,  in  which  our  boundary  value  problem  for  the  Tricomi 
equation  is  formulated,  is  much  simpler.   In  this  space  the 
complex  characteristic  coordinates  are  given  by 

C  =  X  +  iY, 
and 

n  =  X  -  iY  . 

The  sonic  surface  corresponds  to  the  complex  plane  Y  =  0. 

For  analytic  continuation  around  Y  =  0,  a  path  that  connects 

the  points  (X,Y)  =  (0,1)  and  (X,Y)  =  (0,-1)  can  simply  be 
defined  by 

X  =  0  ,      Y  =  e^^  ,      0  <  e  <  TT  . 

The  initial  characteristic  planes  through  the  point 
(0,1)  are  given  by  C  =  i  and  n  =  -i.   The  corresponding 
5-path  and  n-path  are  respectively 

^   =    iY   =   ee         ,  0<_6<_iT, 

and 

T-rri    /O      •;  ft 

0      <      6      <      TT      . 
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In  Figure  1  we  present  these  paths.   They  make  it  possible 
to  extend  solutions  of  (1.2)  from  the  upper  unit  semicircle 
in  the  X,Y  plane  through  the  complex  domain  into  the  lower 
unit  semicircle. 
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3.    GREEN'S  AND  NEUMANN'S  FUNCTIONS  AND  THE  POISSON  FORMULAS 

In  this  chapter  we  give  closed  formulas  for  the  solutions 
of  the  Dirichlet  and  mixed  Dirichlet  and  Neumann  problems 
in  the  upper  unit  semicircle   for  the  more  general  equation 


(3.1)  "XX  -^  '^YY  +  I  '^Y 


with  values  of  k  to  be  specified.   This  is  the  Euler-Poisson- 
Darboux  equation  and  is  also  known  in  the  literature  as  the 
generalized  axially  symmetric  potential  equation  [27]. 

For  k  a  positive  integer  (3.1)  is  the  equation  of  a 
potential  with  axial  symmetry  in  space  of  dimension  n  =  k+2. 
Weinstein  has  considered  equation  (3.1)  as  the  equation 
for  the  axially  symmetric  potential  in  a  fictitious  space 
of  dimension  k+2,  with  k  positive  and  real  [25-28] .   The 
parameter  k  will  be  referred  to  as  the  index  of  the 
equation. 

A  fundamental  solution  valid  for  values  of  k  >  0  is 
given  in  Section  3.1.  This  fundamental  solution   is  the 
Neumann  function  for  the  upper  half-plane.   In  Section  3.2 
we  obtain  the  Green's  function  for  the  upper  unit  semicircle 
when  0  <  k  <  1.   In  Section  3 ,3  we  use  this  Green's  function 
to  obtain  a  Poisson  formula  for  the  unit  semicircle; 
that  is,  we  give  a  representation  of  the  solution  with 
vanishing  boundary  values  on  Y  -  0  in  terms  of  its  values 
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on  the  boundary  of  the  semicircle.  In  Section  3.4  we  obtain 
for  k  >  0  the  representation  of  the  solution  when  its 
values  are  prescribed  on  the  boundary  of  the  unit  semicircle 
and  the  normal  derivative  vanishes  on  the  X  axis.   We  obtain 
these  two  representations  independently.   However,  a 
"correspondence  principle"  [26]   that  will  be  used  in 
Section  3.2  relates  solutions  of  (3.1)  with  different 
indices  and  makes  it  possible  to  deduce  the  solution  of 
the  Dirichlet  problem  from  that  of  the  Neumann  problem. 


3 . 1   The  Fundamental  Solution 

A  fundamental  solution  for  (3.1)  with  index  k  >  0  is 
given  by 


(3.2)   S(X,Y;C,Ti)  =  n 


u  •  k-1 

k  sin    a  da 


[(X-U^+  (Y-n)^+  4Yn  sin^  I  ]''/2 


where  the  singularity  is  at  the  point  (C/H).   This  is  a  general- 
ization of  the  Gauss  formula  for  the  gravitational 
potential  of  a  circular  wire  in  a  three-dimensional  space, 
to  which  it  reduces  if  we  make  k  =  1  ([10],  p.  148). 
Similar  formulas  have  been  given  by  Weinstein  [25-28] , 
although  there  is  some  confusion  in  the  literature. 
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Fundamental  solutions  have  been  known  for  the  Tricomi  equa- 
tion since  the  time  of  his  original  paper  [23] . 

It  can  be  checked  by  direct  calculation  that  (3.2)  is 
a  solution  of  equation  (3.1).   Moreover,  as  a  function  of 
the  parameters  C  and  n  ,  (3.2)  is  a  solution  of  the  adjoint 
equation 


,3.3)         vj^  .v^^  -^v^  .^v=0  . 


We  introduce  the  notation 


pj  =  (x-C)^  +  (Y-n)^  ,   qI  =    (X-?)^  +  (Y+n)^  , 


(3.4) 


2     ^l 


q  = 


^1 

with  this  notation  (3.2)  can  be  written  in  the  form  [25] 


n  k/2  f      k-1      2  ('^-2)/2r   2       . -k/2 

(3.5)  S(X,Y;5,n)  =  (§)^     t*^    (1-t^)         [^^  ■"  ^  J     '^^ 

This  is  a  hypergeometric  integral.   It  can  be  expressed  as 

~k-l   k   /-p  ,k,  ^  2 

(3.6)  S(X,Y;C,ri)   =   ^ —rjk) F||,|;k;l-q' 

by  use  of  the  relationship  [18] 
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1 

(3.7)   F(a,b;c;z)  =  r  (brr)c-b)  j  ^^~^  d-t)""'^"^  (1-zt)"^  dt, 

0 


valid  for  Re  c  >  Re  b  >  0  ,  |arg(l-z)  |  <  tt  . 

The  h7pergeometric  differential  equation  has  three  singular 
points,  and  at  each  of  those  points  one  of  the  two  linearly 
independent  solutions  has  a  logarithmic  singularity.  Using 
known  relations  for  the  hypergeometric  function  [17],  we 
can   rewrite  formula  (3.6)  as 


(3.8)  S(X,Y;5,n)  =  ^-^   |F(|,|;l;q^]  log  i  +  G(X,Y;?,n)}  , 


where  G  is  a  regular  function.   The  coefficient  of  the 
logarithmic  term,  namely, 

(3.9)        A(X,Y;C,n)  =  ^-g-  F[|,|;l;q^ 

is  the  Riemann  function  of  equation  (3.1),  which  shows 
that  (3.8)  has  the  correct  singularity  for  a  fundamental 
solution. 


17 


3. 2   The  Green's  Function  for  the  Unit  Semicircle 

Let  z(k)  be  a  solution   of  equation  (3.1),  where 
temporarily  we  let  the  index  k  be  either  positive  or 
negative.   The  following  "correspondence   principle" 
has  been  introduced  by  Darboux  and  extensively  used  by 
Weinstein  and  Diaz  [7,26]: 


(3.10)  z(k)  =  Y-"-"^  z(2-k)  . 


The  proof  is  done  by  direct  verification.   For  0  <  k  <  2 
it  relates  solutions  with  positive  index.   Using  this 
principle  we  have  from  (3.2)  a  new  fundamental  solution  [28; 


7T  ^  , 

(3.11)   U(X,Y;?,n)  =  nY^-'^  [  2 ^"   "  ^"^ 


Q  [(X-a''+(Y-n)'+4Yn  sin^  I  ]  ^2  '^^/^ 


The  results  of  Section  3.1  applied  to  the  solution  of  index 
(2-k)  show  that  (3.11)  has  the  proper   singularity.  For 
0  <  k  <  1   the  expression  (3.11)  is  a  Green's  function 
relative  to  the  upper  half -plane,  since    it  vanishes  on 
the  X  axis. 

We  use  the  method  of  inversion  to  obtain  the  Green's 
function  for  the  unit  semicircle.   Our  idea  is  to  extend 
Kelvin's  transformation  to  equation  (3.1).   More  precisely, 
it  can  be   proved  by  direct  calculation  that  if  u(X,Y) 
is  a  solution  of  (3.1),  then  so  also  is 
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u'  (X,Y)  =  r  '^  u(X'  ,Y' 


where 

X'  =  -^   ,        Y'  =  ~   ,  r^  =  X^  +  Y^ 

r  r 


Applying  Kelvin's  transformation  to  (3.11),  we  obtain  a 
new  solution  of  (3.1)  with  respect  to  the  variables  X  and 
Y,  namely, 

(3.12)     U'(X,Y;^,n) 


1-k 

1-k  k-2  \   sin    a  dg 


We  recall  that  (3.11)  is  a  solution  of  the  adjoint 

equation  (3.3)  as  a  function  of  the  variables  ?  and   n . 

Using  the  fact  that  if  u(X,Y)  is  a  solution  of  (3.1), 
we  find  that 


(3.13)  v(X,Y)  =  Y^  u(X,Y) 


is  a  solution   of  (3.3).   Applying  Kelvin's  transformation, 
we  obtain  the  new  solution 

(3.14)     V(X,Y;C,n) 


1  -k 

1-k    k-2   '  -•   -L  K 

nY    p 


sm g  da 

i  [(X-5')^+(Y-n')^+4Yn'  sin^  ^ 


1 (2-k)/2 
2  ^ 


19 


of  (3.1),  where 


P  P 


Formula  (3.14)  defines  a  fundamental  solution  of  (3.1) 
with  singularity  at  the  point  (?',n'),  which  is  the 
inverse  image  of  (C^n)  with  respect  to  the  unit  circle, 
It  can  be  verified  that  (3.11)  and  (3.14)  take  on 
the  same  values  on  the  boundary  of  the  unit  semicircle, 
For  that  we  have  just  to  notice  that  for  (X,Y)  on  the 
circumference  of  radius  1,  the  relation 


^=    R 

P 


holds,  where 


R^  =  (X-O^  +  (Y-n)^  ,     R^  =  (X-C)^  +  (Y-n')^ 


Subtracting  (3.14)  from  (3.11),  we  obtain  the  Green's 
function  for  the  unit  semicircle.   Let  us  introduce  the 
abridged  notation 


(3.15)   I  (X,Y;5,n)  =  I  2 2 2"  "^  p/2    * 

P  i  [(X-C)  +(Y-n)  +4Yn  sin^  I  ]P/^ 

Then  the  Green's  function  for  (3.1)  in  the  unit  semicircle, 
with  index  0  <  k  <  1,  is 


20 


(3.16)   G(X,Y;5,n)  =  HY^  ^  [1^_^{X  ,Y  ;  g  ,r]) -p^    ^1^_^(X  ,Y  ;  Z'  ,r\' ) 


In  the  next  section  we  shall  insert  this  result  into 
Green's  formula  to  obtain  a  Poisson  formula  for  the  unit 
semicircle. 


3. 3   The  Poisson  Formula  for  the  Unit  Semicircle 

The  Green's  function  obtained  in  Section  3.2  can  be  used 
to  solve  the  Dirichlet  problem  in  the  unit  semicircle. 
We  want  to  obtain  the  solution  to  (3.1)  which  takes  given 
continuous  values  on  the  upper  boundary  9D  of  the  unit 
semicircle  D  and  which  vanishes  at  Y  =  0 .   We  start  from 
the  representation  formula  [10] 


(3.17)   u(X,Y)  =  2¥   j  ^(^'^)  Iv  ^M^^'"^'^'^^  ^^ 
9D 


9/9 V  stands  for  the  inner  normal  derivative , and  s  is  the 
arc  length. 

For  G„  we  have  the  interchange  rule 


Gj^(C,ii;X,Y)  =  Gj^(X,Y;C,n) 
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where  G   is  the  Green's  function  (3.16)  of  (3.1).  Formula 
(3.17)  can  be  applied  to  any  domain  such  that  the  coeffi- 
cients of  equation  (3.1)  are  regular  up  to  the  boundary. 
In  particular,  we  can  apply  it  to  a  domain  D   obtained  by 
deleting  the  strip  y  <  e  from  the  unit  semicircle.  By 
continuity  we  can  then  extend  the  formula  to  the  case  in 
which  D  is  the  unit  semicircle  itself.   We  can  then  write 

(3.17)  as 

(3.18)  u(X,Y)  =    -    ~         I  f  (i)  ~   ■(y-^"'^p  sin  (t)[I,  ,  (X,Y;C,ri) 


0 
P^"^  l2_j^(X,Y;?',n')]l  d*  , 


where  f  ((}>)  represents  the  values  of  the  solution  u  on  the 
boundary  of  the  semicircle,  and  where  use  has  been  made  of 
the  fact  that  the  solution  vanishes  at  Y  =  0 .   In  this 
context  we  have  introduced  the  notation 

^  =  p  cos  (J)  ,    n  =  P  sin  cj)  . 

The  kernel  in  (3.18)  can  be  computed  by  a  straightforward 
calculation.   We  obtain,  using  the  notation  (3.15), 


(3.19)   1^  =  (2-k)Y-'  "^    sin  <^[1^_^{X  ,Y  ;  ^  ,n)    +    2B1^_^{X  ,Y  j  ^  ,t])  ] 


at  p  =  1 ,  where 


2  ct 
X  cos  4)  +  Y  sin  ())  -  1  -  2Y  sin  (J)  sin   y 
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Manipulating  formula  (3.19),  we  obtain 


(3.20)      11=  (2-k)  Y^  ^  sin  <i>(r^-l)    I^_j^  (X ,  Y;  Cl ) 


Therefore, the  solution  (3.18)  can  be  written  as 


IT 


0 


(3.21)   u(X,Y)  =  ^  I  f  (0)  (l-r^)Y-^  ^ 


J  (   sin g  da 1 

•  =^"n  J  ,(x-5)^(v-n)2MYnsin2|l<^-'^'/2/ 


when 


E,   =   cos  4)  ,    n  =  sin  (j)  . 

The  hypergeometric  integral  under  the  bracket  can  be 
transformed  into  a  hypergeometric  function  using  results 
similar  to  those  presented  in  Section  3.1.   We  obtain  for 
this  integral  the  formula 

kx^2 


"^^^^^^r^t-i--!'--^'--^)- 


Making  use  of  the  quadratic  transformation 


P„.  |,  ,.  |,  .-.,.,  =  a-z,-<l->^/^'  rfl-  |,  -  |;|  -  |,j4| 


for  the  hypergeometric  function  ([18],  p.  51),  we  obtain 
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from  (3.21)   after  further  transformations,  the  representa- 
tion 


(3.22)  u(X,Y)  =  CY^-^     f(cD)[l-(X^Y^ 
J   n+y  +Y  -9y  <-n 


sin  (j) 

(4-k)/2 


where 


i7ii    k   3    k   3    k        4Y   sin  4)     ],, 

■  ^P  "  4 '  2  "  4'  2  "  2'  2 2 ^ 2^^*' 

*■     4   z    4   2    z   (1+X^+Y^-2X  cos  (i))^J 


(r(i  -  \)]^ 


2\  ^(2-k) 


This  formula  generalizes  a  result  of  Germain  and  Bader  [8] 
that  was  given  for  the  case  k  =  ^  ,  and  is  the  representa- 
tion we  have  been  seeking.   It  solves  the  Dirichlet  problem 
for  the  Euler-Poisson-Darboux  equation  with  index  0  <  k  <  1 
in  the  unit  semicircle  and  with  vanishing  values  of  the 
solution  at  Y  =  0. 
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3 .4   The  Mixed  Dirichlet-Neumann  Problem 

The  second  problem  we  are  interested  in  is  the  mixed 
boundary  value  problem  in  which  the  values  of  the  solution 
are  given  on  the  upper  boundary  of  the  semicircle  and  the 
normal  derivative  vanishes  at  Y  =  0 .   We  ask  that  the  solu- 
tion be  continuous  up  to  the  boundary.   We  have  already 
found  in  (3.2)  a  Neumann's  function  for  the  upper  half -plane, 
valid  for  k  >  0.   It  can  easily  be  checked  that  the  normal 
derivative  of  this  function  vanishes  at  Y  =  0 . 

In  order  to  obtain  a  representation  formula  for  the 
solution  we  have  to  devise  a  formula  different  from  (3.17), 
since  we  are  dealing  with  a  Neumann  condition  at  Y  =  0 . 
The  difficulty  is  that  we  have  an  equation  which  is  not 
self-adjoint,  so  we  have  to  find  a  function  that  has  simul- 
taneously the  properties  of  a  Neumann's  function  for  our 
equation  and  of  a  Green's  function  for  the  adjoint  equation. 

We  start  with  the  identity  [10] 

(3.23)  \^i^        (vL[u]-uM[v])  dX  dY  +  I  [v  ^  -u  |^  +(a  |^  +b  |f)uv]ds 
D  9D  ^  n 


where 


L[u]  =  u^^   +   Uyy  +  a(X,Y)Uj^  +  b(X,Y)uY  +  c(X,Y)   =   0 


and  where   M[v]  =0   is  the  adjoint  equation, 
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This  formula  is  valid  if  the  coefficients  a,  b,  and  c  are 
regular  up  to  the  boundary  of  the  domain  D.  In  our  case 
we  have 


a=0,    b  =  |,    c  =  0 


Formula  (3.2)  can  be  written  as 


^L  "  "^^    Jj^(X,Y;C,n)  , 


where  we  use  the  notation 


■^  k-1 

sin    a  da 


J  (X,Y;C,n)  -  I    -  ^5 T ^ j^      . 

P  I    [(X-^)^(Y-n)'  +  4Yn  sin^  |]P/2 


Then, 


N  =  Y*"  J^(X,Y;C,n) 


is  a  fundamental  solution  for  the  adjoint  equation 


^f^^  =  -XX  ^  -YY  -  I  -Y  ^  ^  -  =  ° 


By  use  of  the  Kelvin    transformation  we  obtain,  as 
in  Section  3.3,  the  Green's  function 


(3.24)   Gj^(X,Y;C,n)  -  y'^  [  Jj^  (X  ,  Y;  5  ,n  )  "  P  '^J^  (X  ,  Y  ;  ^  '  ,n  '  )  ] 
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where 


C        ,     n         2..   ^2     2 


The  function  G.,  vanishes  on  the  whole  boundary  of  the  semi- 

M 

circle,  and  the  expression  in  square  brackets  in  (3.24), 
is  a  solution  of  L[u]  -  0  which  has  vanishing  derivative 
with  respect  to  Y  at  Y  =  0 . 

We  apply  the  identity  (3.23)  to  a  solution  u(X,Y)  of 
(3.1)  and  to  v  =  G., .   We  take  as  domain  of  integration  D   / 

M  £ 

the  unit  semicircle  with  a  strip  0  <_  Y  <^  e  deleted.  Since 
G   behaves  like  a  Green's  function  on  the  upper  boundary 
r  of  the  semicircle,  we  obtain 

2.u(e,n)  =  I  u(X,Y)  I-  G^  ds  -  J  [g^  1^  -  u  g^  +  I  ugJ  dX 


where  we  note  by  y  the  segment  Y  =  e  of  D  .   The  second 
integral  gives,  after  a  cancellation  and  putting 
Jj^  -  Jj^(X,Y;C'  ,n')  , 


3G^ 

dX 


1  K  If  -  "  av^  *  I  "=m] 


Y  *- 
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If  we  now  let  e  ^  0,  the  whole  integral  over  y  vanishes 
once  we  assume  8u/8Y  bounded.   Thus  we  obtain  the  formula 


u(C,n)  =  -  2¥  f  ^^*)  1^  Y^[Jj^(X,Y;5,n)  -  P  ''jj^(X,Y;^',n')]  d4). 


0 
or,  equivalently , 


u(X,Y)  =  -  2?  I  ^(*)  J^   n'^[Jj^(X,Y;C,n)  -  r  "^Jj^  (X  '  ,  Y  '  ;  C,  H )  ]  d(}). 


where 


-^  ,    Y'  =  -|-  ,    and     r^  =  X^  +  Y^ 


Following  the  same  steps  as  in  Section  3.3   we  are  led 
to  the  formula 

•T  ^2      k-1 

^v  w\     k  f  .  k,  c,..(    f       (1-r  )  sin     a  da           ,,, 

u(X,Y)  =27    sm  (I)  f((}))^    2 2 2  a  (k+2)/2^^'^ 

^^  J)  ^  J,  [(X-U%(Y-n)^+4Yn  sin^  a^(k+2)/^/ 


where 


E,   =   cos  4)  ,    n  =  sin  4)  . 

This  formula  can  be  written,  by  using  similar  relations  for 
the  hypergeometric  function,  as 
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f((j))  [l-(x^ 
2   2 
J   (1+X  +Y  -2X 


(3.25,  u(X,V,  =  c  I   --  -  -'^^'''  t^;2?/2 

Q   V-.-   .  -   —  COS  (J))  '     '  ' 


with 


,-k    1   k     k    1,  4Y   sm  4) -, 

^*-4    2'  4   -•-'  2    2'  2 2 o"-'^*' 

^    ^   ^      ^    ^   (1+X'^+Y^-2X  cos  ^)^ 


,k-2  (r(^))^ 
7T    r{k) 


The  result  is  valid  for  k  >  0.   It  is  the  desired  represen- 
tation of  the  solution  of  the  mixed  Dirichlet-Neumann  problem 
in  the  unit  semicircle,  with  vanishing  normal  derivative  on 
the  X  axis,  in  terms  of  the  values   of  the  solution  on  the 
boundary  of  the  semicircle. 

The  correspondence  principle  (3.10)  can  be  applied  to 

the  solution  (3.25)  with  index   k'  =  2-k   and  boundary  values 

k-1 

sin    (})f((t)),   0<k<l.   We  obtain  from  it  the  solution 

(3.22)  of  the  Dirichlet  problem  with  boundary  values  f  ((})), 
as  was  mentioned  earlier. 
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4.    THE  REFLECTION  PRINCIPLE 

We  establish  in  this  chapter  reflection    laws  for 
solutions  of  the  general  Euler-Poisson-Darboux  equation 


(4.1)  u^^  +  u^Y  +  I  u^  =  0  ,  0  <  k  <  1 


which  make   it  possible  to  extend  them  from  the  upper  to  the 
lower  unit  semicircle.   They  will  be  used  in  Chapter  5  to 
prove  our  main  theorem  about  the  transonic  boundary  value 
problem.   However,  they  have  an  interest  of  their  own. 
General  theorems  about  reflection  laws  for  solutions  of 
elliptic  partial  differential  equations  with  analytic  condi- 
tions on  the  boundary  are  well  known  (Garabedian  [11],  Lewy 
[16]).   The  special  feature  of  our  case  is  that  the  boundary 
over  which  the  reflection  is  going  to  be  made  is  a  singular 
line  for  the  equation. 

It  is  known  (Hyman  [12])  that  a  solution   of  (4.1) 
which  is  analytic  in  some  neighborhood   of  a  segment  of  the 
X  axis  is  necessarily  even  in  Y  and  is  uniquely  determined 
by  its  values  on  this  segment.   This  is  a  local  theorem 
based  on  power  series  expansion  of   the  analytic  solution, 
valid  for  k  different  from   0   or  a  negative  integer.  The 
result  for  k  =  1  is  classic  [13].   However,  this  theorem  does 
not  apply  to  solutions  with  a  branch  singularity  of  the  type 
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(4.2)  u(X,Y)  =  Y-"-  ^  f{X,Y^) 


Our  representation  formula  (3.22)  is  an  example  of  this  type 
of  solution.  We  shall  derive  reflection  laws  that  are  valid 
in  the  large  and  permit  one  to  reflect  solutions  of  the  type 
(4.2)  for  which  the  X  axis  is  a  branch  line.  We  use  analy- 
tic continuation  of  the  solutions  into  the  complex  domain 
of  the  variables  X  and  Y  to  do  this. 

Let  the  variables  X  and  Y  be  extended  to  complex  values 

Y  =  Y^  +  iY2  =  |Y|e^®  . 

We  define  the  path  y   by 

Y  =  I  (X,Y:  X  =  0,  Y  =  re^^,  0  <  6  £  it  i, 

where  0  <  r  <  1.   Consider  the  analytic  continuation  of  the 
solutions  (3.22)  and  (3.25)  along  this  path  y.      The  explicit 
form  of  the  solutions  guarantees  that  the  continuation  can 
be  performed.   Two  points  A  and  B,  with  A  in  the  upper  unit 
semicircle  and  B  in  the  lower,  can  be  connected  by  a  path  of 
the  form 


(4.3) 


where  y,       is  a  path  contained  in  the  upper  semicircle  which 
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connects  the  points  A  and  (0,r),  and  -^ ^    is  the  corresponding 
path  which  connects  B  and  (0,-r)  in  the  lower  semicircle. 
By  the  monodromy  theorem,  any  other  path  which  does  not  wind 
around  the  X  axis  will  give  the   same  result. 

We  obtain  in  this  fashion  the  following  theorem:  If 
u(X,Y)  is  a  solution  to  the  Dirichlet  problem  for  equation 
(4.1)  in  the  upper  unit  semicircle  with  zero  boundary  values 
u  =  0  at  Y  =  0,  then  it  can  be  extended  to  the  lower  semicircle 
by  the  reflection  law. 


(4.4)      u(X,-Y)  =  e^-"-"^^^^  u(X,Y)  ,     Y  >  0 


where  analytic  continuation  has  been  made   along  the  path 
(4.3).   Similarly,  the  solution  of  the  mixed  Dirichlet  and 
Neumann  problem  in  the  same  domain,  with  vanishing  normal 
derivative  ^y   ^   0  at  Y  =  0,  can  be  extended  to  the  lower 
semicircle  according  to  the  law 

(4.5)  u(X,-Y)  =  u(X,Y)  ,      Y  >  0  . 

Clearly,  the  reflection  rules  depend  only  on  the  form  of  the 
boundary  condition  at  Y  ==  0  and  are  valid  in  more  general 
domains. 
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5.    SOLUTION   OF  THE  BOUNDARY  VALUE  PROBLEM 

We  can  now  state  and  prove  our  main  theorem.   Consider 
the  Tricomi  equation 


(5.1)  yu    +  u 

■^  XX     yy 


and  the  transformed  equation 


(^•2)  ^XX  -^  ^YY  -^  sf-^  =  °  ' 


in  canonical  coordinates  X  and  Y.   We  introduce  the  class  of 
analytic  solutions  of  equation  (5.1)  which  are  single-valued 
in  a  neighborhood   of  the  line  y  =  0.   Let  u(X,Y)   be  a 
real  solution  of  (5.2)  in  the  upper  unit  semicircle 


v  =  f 


D^  =  -1  (X,Y)  :  X^  +  Y^  <  1  ,  Y  >  0 


belonging  to  this  class.   Then  in  the  complex  domain  of  the 
extended  variables  X  and  Y  the  solution  u(X,Y)  can  be  continued 
analytically  along  a  path,  as  described  in  (4.3),  into  the 
lower  unit  semicircle. 


D_  -  |(X,Y):  X^  +  Y^  <  1  ,  Y  <  ol 


We  shall  prove  that  the  boundary  value  problem 
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Re  ju(X,Y)^   =  F(e)  ,    X  =  cos  6  ,   Y  =  sin  6  ,   0  <  9  <  2tt, 

where  F(9)  is  periodic  and  continuous,  is  well  posed  within 
the  specified  class  of  solutions  u(X,Y). 

We  give  first  an  explanation  of  the  choice  of  the  class 
of  functions  that  we  have  made.   If  u(x,y)  is  single  valued 
and  analytic,  it  admits  a  double  power  series  expansion  in  a 
neighborhood   of  the  origin  of  the  form 


/    \    V      n  m 
u(x,y)  =  Z  a^^  X   y   . 


Therefore. we  can  write 


~      3      ~      3      2~      3 

u(x,y)  =  V, (x,y  )  +  yv2(x,y  )  +  y  v^{x,y    )  , 


u(X,Y)  =  v^(X,Y^)  +  Y^/^  V2(X,Y^)  +  Y^^^    v^{X,Y^)     , 

where  v,  ,  v„  and  v^  are  regular. 

If  we  let  the  variable  Y  wind  around  the  X-axis,  through 
an  angle  of  2-n   or  4tt,  we  obtain  new  solutions 


u'(X,Y)  =  v^(X,y2)  +  e^^i/3  y2/3  v2(X,y2)  +  e^^i/^Y^/\3  (X, Y^) 
and 


u"(X,Y)  =  v^(X,y2)  -h  e^^^/^Y^/^^^^'^^)  -^  el^^^/\^/\3(X,Y2) 


of  equation   (5.2),  as  the  equation  itself  remains  unchanged. 
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We  apply  now  the  differential  operator  L  defined  by  equation 
(5.2)  to  u ,  u'  and  u" .   This  gives 


L[v^]  +  LlY^^^v^]    +    L[Y^^^v^]    =    0 


L[v^]    +    e^-i/^    L[y2/3v^]    +    e^^'^'    L[y'/\^]    =    0 


L[v^]    -.   e«-i/^   UY^/\^]    +   el^-i/3   L[y4/\3] 


Because  the  determinant  of  this  linear  system  of  equations 

is  different  from  zero,  it  follows  that  each  of  the  terms 

3  3       2       3 

v^  (x,y  ),  yv  {x,y  ),  and  y  v-.(x,y  )   is  a  solution  of  the 

Tricomi  equation  separately.   By  the  uniqueness  theorem  for 

the  Cauchy  problem  for  (5.1),  we  can  conclude  that 


y^v,(x,y^)  E  0  . 


Therefore, u (X,Y)  admits  a  unique  decomposition  of  the  form 


(5.3)      u(X,Y)  =  v^(X,Y^)  +  Y^/^v^(X,Y^) 


We  proceed  now  to  prove  the  main  theorem.   Let  the 
boundary  function  F(e)  be  defined  by 


u(X,Y)  =  h(X)  ;    Y  >  0  ,      X^  +  Y^  =  1 
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and 


ju(X,Y)| 


Re  ■(u(X,Y))-  =  k(X)  ;    Y  <  0  ,     X^  +  Y^  =  1  . 


Let  u,  (X,Y)  be  a  solution  of  the  Dirichlet  problem  in  the 
upper  unit  semicircle  such  that 


2     2 

X+Y=l,      Y>0 


with  f  (X)  continuous  and  vanishing  at  X  =  +  1.   Let  u^(X,Y) 
be  a  solution  of  the  mixed  Dirichlet-Neumann  problem,  in 
the  same  domain,  with 


U2(X,Y)  =  g(x)  ,      X^  +  Y^  =  1  ,      Y  >  0  , 

9u 

y^        (X,0)   =   0   . 


Both  u,  (X,Y)  and  U2(X,Y)  can  be  computed  explicitly  by  means 
of  the  representation  formulas  (3.22)  and  (3.25)  of  Chapter  3. 
By  the  reflection  laws  obtained  in  Chapter  4,  we  have 


u^(X,-Y)  =  e^^^/^  f(X)  ,     Y  >  0  , 
and 
u  (X,-Y)  =  g(X)  ,  Y  >  0  ,     X^  +  Y^  =  1  . 
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If  we  choose  f  and  g  so  that 


f  (X)  +  g(X)  =  h(X) 


f  (X)  cos  4  TT  +  g(X)  =  k(X) 


that  is, if  we  set 


and 


then 


f{X)  =  I  (h(X)  -  k(X)] 


g(X)  =  5^  +  §  k(X) 


(5.4)         u(X,Y)  =  u^(X,Y)  +  U2(X,Y) 


is  a  solution  of  (5.2)  which  satisfies  the   prescribed 
boundary  conditions. 

The  uniqueness  of  the  solution  of  the  boundary  value 
problem  follows  immediately.   The   representation  formulas 
(3.22)  and  (3.25)  show  that   u^(X,Y)  and  U2(X,Y)  have  the 
form 

u^(X,Y)  =  Y^/^v^(X,Y^) 


and 


U2(X,Y)  =   V2(X,Y^) 


respectively.   The  decomposition  (5.4)  is  therefore  of  the 
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form  (5.3)  and  is  unique.   The  boundary  condition  on   u  is 
seen  to  impose  Dirichlet  and  Dirichlet-Neumann  data  on  u, 
and  u_  respectively,  which  determine  them  uniquely,  complet- 
ing  the  proof.   The  continuous  dependence  on  the  boundary 
values  is  also  a  consequence  of  the  explicit  representation 
we  have  given  of  the  solution. 

We  remark  that  the  theorem  could  be  extended  to  arbi- 
trary  domains  symmetric  with  respect  to  the  X  axis  by 
performing  reflection  in  the  manner  described  in  this  work, 
but  resorting  to  a  more  general  treatment  of  the  Dirichlet 
and  Dirichlet-Neumann  problem  for  u,  and  u^. 
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